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Abstract. "Cometary" Blue Compact Dwarf Galaxies (iI,C 
BCDs) are characterized by an off-center starburst close to the 
end of their elongated stellar bodies. This rare phenomenon 
may carry some clues on how collective star formation ig- 
nites and propagates in gas-rich low-mass stellar systems. This 
off-center burst may be a fortuitous enhancement of the other- 
wise moderate star-forming activity of a dwarf irregular (dl), 
or may be caused by a set of special properties of such sys- 
tems or their environment. We attempt here a first investi- 
gation of this issue by analysing two prototypical examples 
of cometary dwarf galaxies, the nearby iI,C BCDs Markar- 
ian 59 and Markarian 7 1 , both containing an extraordinarily lu- 
minous Hll region in the outskirts of a dl-like host. Using 
deep ground-based spectrophotometric data [}| and HST im- 
ages f], we study the physical state of the starburst regions and 
the structural properties of the underlying irregular galaxies. 
We find that the average metallicities show small scatter in the 
vicinity of the star-forming regions and along the major axis 
of Mkn 59 which suggests that mixing of heavy elements must 
have been efficient on scales of several kpc. The azimuthally 
averaged radial intensity distributions of the underlying host 
galaxies in either iI,C BCD can be approximated by an expo- 
nential law with a central surface brightness and scale length 
that is intermediate between typical iE/nE BCDs and dwarf 
irregulars. Spectral population synthesis models in combina- 
tion with colour magnitude diagrams and colour profiles yield 
a most probable formation age of ~ 2 Gyr for the low surface 
brightness (LSB) host galaxies in both iI,C BCDs, with upper 
age limits of ~ 4 Gyr for Mkn 59 and ~ 3 Gyr for Mkn 71, 
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i.e. significantly lower than the typical age of several Gyr de- 
rived for the LSB component of iE/nE BCDs. These findings 
raise the question whether iI,C systems form a distinct physical 
class within BCDs with respect to the age and structural prop- 
erties of their hosts, or whether they represent an evolutionary 
stage connecting young iO BCDs and "classical" iE/nE BCDs. 
In spite of the scarcity of available data, a review of the prop- 
erties of analogous objects studied in the local universe and at 
medium redshifts provides some support for this evolutionary 
hypothesis. 
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ies: evolution — galaxies: compact — galaxies: starburst — 
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1. Introduction 

Despite the effort that has been devoted to the investigation of 
Blue Compact Dwarf Galaxies (BCDs), the origin of their re- 
current starburst activity as well as its impact on their spec- 
trophotometric and dynamical evolution is still poorly under- 
stood. A thorough investigation of these processes is required 
before questions pertaining to evolutionary connections among 
dwarf irregulars (dls), dw arf ellip ticals (dEs) and BCDs (Thuan 

1996b; 



1985L Davi es & Phillipps |1988| Papa deros et al. 1996b; here 
after ^96bj Pa tterson & T huan |1996] , Salzer & Norton |l998 



Marlowe et al. 1997 , 1999 ) can be answered. Deep imaging of 
the low surface brightness (LSB) component of BCDs, discov- 
ered first by Loose & Thuan (1985; hereafter |LT85 ), disclosed 
an evolved stellar population underlying the regions of active 
star formation (SF) (Loose & Thuan 1986, Kunth et al 



1988 



Papaderos et al. 1996a; hereafter ^96a ). Therefore, the initial 
hypothesis that BCDs are "extragalactic H II regions" lacking 



an older stellar population (Sargent & Searle 1970) had to be 
dismissed for the majority of these systems. 

An exception to this finding is made by a tiny fraction (^ 
1% ) of young galaxy candidates wh ich re cently have been 
identi fied among BC Ds (Th uan et al. 1997 , Papaderos et al. 



1998| , Thuan & Izotov |1999J). There is increasing observational 
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evidence that such systems, invariably found among the most 
metal-deficient (Z < 1/20 Z Q , Thuan & Izotov |1999| Izotov 



Table 1. Ground-based imaging 



& Thuan 1999) members of the BCD class, have started form- 
ing the bulk of their stellar component less than 100 Myr ago. 
This qualifies them as local counterparts of primeval galaxies, 
i.e. the first galactic building blocks which are thought to have 
formed at high redshifts. 

It seems meaningful to search among the variety of dwarf 
galaxies for possible successors of young galaxies, i.e. objects 
in a more evolved stage with a substantial fraction of their 
masses still in gaseous form. An age sequence built therefrom 
would allow to study the evolution of nearby dwarf galaxies, 
which may at the same time give clues to early galaxy evolu- 
tion processes at high redshift. This is of great interest for the 
assessment of the starburst-driven evolution of young galax- 
ies, some of which may become the building blocks of larger 
systems and for the understanding of issues like the faint blue 



galaxy excess at intermediate redshifts (Babul & Rees |1992| , 
Koo et al 



1997, Guzman et al. 1998) 



Loose & Thuan (1985) developed a morphological classifi- 
cation scheme for BCDs and found that ~ 90% of their sample 
is made up of iE- and nE-objects. These show, respectively, 
an irregular or nuclear starburst component superposed on a 
smooth elliptical LSB host galaxy with optical colours point- 
ing to ages of several Gyr (Loose & Thuan 1984 P96b, Telles 
& Terlevich 1997, Noeske 1999). In contrast, the few members 



of the rare class of iO BCDs examined so far have been found to 
be young galaxy candidates (Thuan et al. 1997, Papaderos et al. 
1998| , Thuan & Izotov |1999| , Izotov et al. |1999[ ), as they appar- 



ently lack a smooth stellar underlying component that would 
have required previous episodes of SF to be built. 

Only the class of il BCDs, characterized by both irregular 
LSB components and starburst regions, has not yet been studied 
systematically. This paper focuses on the subset of iI,C BCDs 
(LT85) where "C" denotes a cometary appearance caused by a 
bright star-forming complex at one end of an elongated stellar 
galactic body (see Figs. |l] and ||). In the notation of Salzer et 



al. ( |1989[ ) an iI,C BCD is to be included in the class of H II 
hotspot galaxies; these systems seem to systematically differ 
from other subclasses of H II galaxies with respect to the metal- 
licities and excitation properties of their ionized gas. While 
the main bodies of iI,C systems are reminiscent of low sur- 
face brightness dwarf irregulars (dl s) or m agellanic irre gulars 
(Kennicutt et al. |l98C| , Dottori et al. |l994 Wilcots et al. \199$ ) 
with some widespread low-level SF, the bright off-center H II 
regions are typical for BCDs with respect to their sizes, Ha lu- 
minosities and electron temperatures. The detection of Super 
Star Clusters (SSCs) in these spots points to a very r ecent or 
still ongoin g starburst episod e (Kennicutt et al 
etal 



1994,Barthetal. 1994) 



198C, Dottori 



O ptical survey s of BCDs ( |LT85[ Ku nth et al. |1988| , Salzer 
et al. [1989| , p96a[ Marlow e et a l. |1997| , Doublier et al. |l997| ) 
and dls (e.g. Hunter et al. 1998) suggest that in these systems 



SF per unit area correlates with the local surface density of 
the underlying host galaxy. In this respect, the occurrence of a 
starburst with the amplitude observed in iI,C BCDs at one end 
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of an elongated stellar LSB body appears puzzling and may be 
interpreted in at least two ways: 

(i) iI,C BCDs are in fact dls or magellanic irregulars caught 
during a brief stochastic enhancement of their otherwise mod- 
erate star formation rate per unit area. This would imply that 
both the structural properties of their stellar hosts and H I en- 
velopes will be indistinguishable from those of other dls. 

(ii) In iI,C BCDs the conditions necessary for the ignition 
of a starburst are not fulfilled at the center but in the outskirts 
of their LSB components. This may be due to a conspiracy of 
intrinsic and external properties of an evolved BCD (kinemat- 
ics of the gaseous component and Dark Matter distribution or 
external perturbation by a companion or intracluster gas) or al- 
ternatively be a common feature among less evolved systems. 

Here we attempt to test such hypotheses by analysing 
optical data of two nearby iI,C BCDs, Mkn59 (NGC4861) 
and Mkn71 (NGC2366, hosting the Hll region complex 
NGC2363). The distance to Mkn71 was determined to 
£>=3.44Mpc through observations of Cepheids (Tolstoy et al. 
For Mkn 59 different models for infall and small-scale 



1995) 



perturbations in the Virgo Cluster environment yield a distance 
between - 10.7 Mpc (Heckman et al. |1998[) and 17.8 Mpc 



(Tully 1988). The principal results of this work do not depend 



on the exact distance to Mkn 59 (cf. Section 3.1). Throughout 
this paper, we adopt Z?=10.7Mpc (Heckman et al. 1998| ). 

Using broad- and narrow band images and long-slit spectra 
from own ground-based observations and archival HST data, 
we shall investigate the physical state and the chemical com- 
position of the ISM and examine whether the iI,C BCDs under 
study differ from typical iE/nE BCDs with respect to the ages 
and structural properties of their LSB hosts. In Section 2, we 
describe our observations and data analysis. We present our re- 
sults in Section 3, and discuss them in Section 4. In Section 5, 
we summarize our results and conclusions. 



2. Data acquisition and processing 

2.1. Ground-based imaging 
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Fig. 1. (B — I) colour maps with overplotted B band isophotes. left Mkn 59; the isophotes correspond to surface brightness 
lev els of 20, 21, 22 and 22.5-25.5 mag/D" in steps of 0.5 mag. The bright starburst knot (BK), described in detail by Dottori et 
al. ( |1994| > is indicated, right Mkn 71; the isophotes correspond to surface brightness levels of 21, 22, 23, 24, 24.5 and 25 mag/D". 
The giant H II complex NGC 2363, as well as the field used to derive the colour-magnitude diagram of the galaxy's underlying 
stellar population from HST data (CMD 2) are marked. The indices 1-4 along the orientation of the long-slit (centered at the 
axis origin) mark the regions from which the spectra displayed in Figure 8 were extracted. North is up and east to the left. 



2.1.1. Observations and data reduction 

Images were taken on March 7th - 10th 1997 at the 2.2m 
telescope of the German-Spanish Astronomical Center, Calar 
Alto, Spain. We used the Calar Alto Faint Object Spectrograph 
(CAFOS), equipped with a 2048x2048 pixel SITe#ld CCD. 
The focal ratio of f/4.4 in the RC focus and the pixel size of 
24 /im/pixel yield an instrumental scale of 0.53"/pixel and a 
usable field of view of ~ 15'. At a gain ratio of 2.3 e~ ADU^ 1 
the read-out noise was < 3 counts (rms). Column (1) of the ob- 
serving log (Table [[]) contains the equatorial coordinates of the 
targets, cols. (2) and (3) give the filters and exposure times, re- 
spectively. Column (4) lists the night of the observing run, start- 
ing with night 1 from March 7th to March 8th 1997. The seeing 
and the mean sky surface brightness during each exposure are 
given in cols. (5) and (6). During each night, dark-, bias-, and 
flat-field exposures were taken, and th e photometric standard 
field NGC 2419 (Christian et al. |1985| ) was observed at differ- 
ent zenith angles. Using the ESO MIDAS^j software package, 
standard reduction and calibration steps were applied to the raw 
images. During nights 3 and 4 the conditions were photomet- 
ric, resulting in calibration errors well below 0.05 mag. A poor 
atmospheric transparency during night 1 led to strong airmass- 
dependent terms for the B band exposure of Mkn 7 1 . Therefore 
the calibration was accomplished by comparing aperture pho- 
tometry of bright isolated point sources on HST- and ground- 
based images. The photometric uncertainty of the latter B band 
exposure was found to be 0. 1 mag. All exposures for one ob- 



ject were aligned to each other using positions of point sources; 
the transformed images were generated by re-sampling the 
original image using a flux-conserving routine. Colour maps 
were derived after the resolution of different frames had been 
equalized by convolving the exposure with the better resolution 
with a normalized gaussian distribution of adequate width. 

The continuum level in the Ha images was inferred by scal- 
ing the R band images by an empirically determined factor C 
so that the fluxes of field stars were matched between the raw 
Ha- and the scaled R exposures. In turn, the scaled R band 
exposures were subtracted from the raw Ha images, giving the 
emission flux F om at each pixel. The latter approach is strictly 
correct only when the Ha emission line is contributing a minor 
fraction of the line-of-sight R band flux. In regions where Ha 
emission makes a substantial fraction of the photons received 
in the R band the empirical scaling may lead to a slight over- 
estimation of the continuum flux. In these cases the correct Ha 
emission FSL. is obtained from F em as: 



F c 



l-C-T R (RaY 



-F„, 



(1) 



4 Munich Image Data Analysis System, provided by the European 
Southern Observatory (ESO). 



where C is the empirical factor described above and Tfj(Ha) 
the mean transmission of the R band filter at the Ha wave- 
length. The resulting narrow band frames were calibrated us- 
ing continuum-subtracted Ha exposures of the planetary neb- 
ula NGC 2392 for which aperture flux measurements are given 
by Kaler ( |1983| ). 

To derive Ha equivalent width maps we subtracted first the 
corrected line emission frame F£ m (Eqjj]) from the raw Ha im- 
age to obtain a continuum image F cont , which was then nor- 
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Fig. 2. Logarithmic representa- 
tion of the Ha equivalent width 
(EW(Ra)[A}) map of the Hll 
complex NGC 2363 within Mkn71, 
computed from HST/WFPC2 images 
as described in Section 2.1.1. B band 
contours at 18-20 mag/D" mark the 
positions of the young star clusters 
NGC 2363 1 and II within NGC 2363. 
The arc-like structure, extending out 
to ~ 0.5 kpc northeast of NGC 2363 II, 
probably marks a supergiant shel l 
(SGS) (Hunter & Gallagher |1997[ ). 
The expanding superbubble (B) and 
the probable outflow chimney (Ch, 
delimited by the pair of dashed lines) 
around NGC 2363 1 (Roy et al. |1991|) 



are indicated. The spatial distribution 
of the young stellar background in 
NGC 2363 1 and NGC 2363 II is re- 



viewed in Section 4. 1 . North is up and 
east to the left. 



malized to the continuum flux per lA wavelength interval by 
dividing it by the effective width of the Ha filter. The emission 
line image F£ m was divided by the transmission of the Ha filter 
at the Ha wavelength. The resulting frame, containing the total 
Ha emission flux, was divided by the 1A- normalized contin- 
uum flux frame to obtain the Ha equivalent width image. 

The colour- and Ha maps were c orrect ed for interstellar 
extinction following Savage & Mathis ( 1979 ) and using the ex- 
tinction coefficients C(H/3) given in Section 23, which trans- 
late into E{B-V)=0.06 mag and E(B-V)=0.08 mag for Mkn 59 
and Mkn 71, respectively. 

2.1.2. Surface photometry 

Surface brightness profiles (SBPs) are the result of a transfor- 
mation of a galaxy's 2-dimensional flux pattern into a mono- 
tonically decreasing 1-dimensional intensity distribution. For 
a set of methods to compute the photometric radius R*(/j,), i.e. 
the radius of a circle with an area equal to the one enclosed 
by the isophote at the surface brightness level //[mag/D"], see 
e.g. Loose & Thuan (1986) and P96c. The method of sum- 



ming up the area defined by pixels with fluxes exceeding a 



threshold I(p) (method iii; P96a) was considered most appro 



priate for the BCDs analyzed here as both their underlying and 
starburst components are of irregular morphology. Other tech- 
niques, such as isophote integration or ellipse fitting, though 
well suited for more regular BCDs, are hardly applicable to 
the optical images of Mkn 59 and Mkn 71, mainly due to the 
fact that their high surface brightness components split into 



many star-forming regions occupying a substantial fraction of 
the host galaxy's surface. 

For each surface brightness level fi, the number N(fi) of all 
pixels within a polygonal aperture with a count level > I(fi) 
was determined; the corresponding photometric radius is 



7T » 



(2) 



where A px \ is the solid angle per pixel in D". Numerical simu- 
lations showed that for the technique described above Poisson 
photon noise may cause a systematic flattening of SBPs at a 
signal to noise S/N level 10. This frequently observed ar- 
tificial profile flattening could be, however, satisfactorily cor- 
rected down to a (S/N) level of a 1 by using the adaptive 
filtering algorithm implemented in MIDAS (see Richter et al. 
1991 for a description). The uncertainties at each photometric 
radius were calculated taking into account the corresponding 
intensity, the sky noise of each image, and the number of pix- 
els involved in its calculation (cf. ^96a ). Colour profiles were 
computed by subtracting SBPs from each other; the surface 
brightness- and colour profiles displayed in Figure 3 are cor- 
rected for interstellar extinction (cf. Section [2.1.1[ ). 

Unlike stellar systems like globular clusters and giant el- 
lipticals where the mass to light (M/L) ratio may be regarded 
nearly constant over the whole system, the integrated luminos- 
ity of BCDs originates from two distinct stellar populations 
with substantially different ages and M/L ratios. Although the 
LSB component contributes on ave rage one half of the total B 
band luminosity of a BCD ( |P96b| , Salzer & Norton 1998), it 
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Fig. 3. left Surface brightness profiles of Mkn 59 (top) and Mkn 71 (bottom) in B, R and /. Small circles show the surface 
brightness distribution of the residual B band emission in excess of the exponential fit to the LSB component (straight line). The 
calibration uncertainty <r ca i effecting a vertical shift of the entire profile is displayed 10 times enlarged, right Radially averaged 
colour profiles with 2 times enlarged calibration error bars. Note that for surface brightness levels fainter than ~ 23.5 B mag 
arcsec -2 the colour index of the underlying LSB galaxy becomes nearly constant. 



contains most of the stellar mass of the system. Therefore, its 
intensity distribution yields information on the inner gravita- 
tional potential within which the starburst ignites and evolves. 
In order to disentangle the light distribution of the older host 
galaxy from the one of the superimposed starburst component, 
we applied a simple 2-component decomposition scheme ad- 
justed interactively to each profile, instead of iteratively fitting 



the full nonlinear 3-component scheme described in P96e . The 



intensity distribution of the LSB component is approximated 
by an exponential fitting law of the form 



R* 



(3) 



or equivalently 

H(R*) = hem + 1-086 



R* 



(4) 



where a denotes its exponential scale length in arcsec and /i 
the surface brightness level in mag/D". Equation (Q) was ad- 
justed to each profile by applying an error weighted linear fit 
to the data, at radii sufficiently large to be free from the con- 
tamination of the starburst light, i.e. where the colour profiles 
become constant and Ha emission vanishes. To check at which 
surface brightness levels gaseous emission becomes negligible, 
we overlaid the Ha emission- and EW(Ha) maps with broad 
band isophotes (cf. Figs, [j] and ||) . The fits, obtained at radii 
R* >76" for Mkn59 and R* > 160" for Mkn71, yield the 
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Table 2. Structural properties of the starburst- and LSB components of Mkn 59 and Mkn 7 1 derived from profile decomposition. 
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extrapolated central surface brightness he,o and the exponen- 
tial scale length a. The radial surface brightness distribution 
of the starburst component can be computed from the residual 
luminosity in excess of the fit Eq. (Q). 

Table [| summarizes the results of the profile decomposi- 
tion. In the same way as described in P96a, cols. (3) and (4) 
give the central surface brightness and exponential scale length 
of the LSB component. Columns (5) and (7) list respectively 
the isophotal radii of the starburst component P25 (practically 
the "plateau radius" used in ^96a ), and of the stellar LSB host, 
£?25- Both radii were determined from extinction-corrected 
SBPs at a surface brightness level of 25 mag/D". Columns (6) 
and (8) contain the apparent magnitudes of the latter compo- 
nents determined within P25 and E25, respectively, and col. (9) 
the total apparent magnitude in each band as derived from in- 
tegration of the corresponding SBP out to the last measured 
point. Column (10) lists the effective radius r e f / and the radius 
r§o, enclosing 80% of the galaxy's total flux. 

2.2. HST images 
2.2.1. Data reduction 

Next we shall discuss colour-magnitude diagrams (CMDs) of 
the Hll complex NGC2363 and the main body of Mkn 71 
derived from Hubble Space Telescope (HST) archival data. 
The following analysis is based on integrations of 2 x 800 and 
1 x700 sec in the F439W filter, 2x800 sec in the F547M filter 
and 1 x 900 plus 1 x 600 sec in the F656N filter taken on January 
8th, 1996 with the WFPC2 (PI: Drissen, Proposal ID 06096). 
All images were reduced through the standard pipeline as de- 
scribed in Holtzman et al. ( 1995a[ ). Exposures taken in the same 
filter were co-added and corrected for cosmic ray events using 
the STSDAS package and IRAF^. We applied a charge transfer 
efficiency correction to the data as described in the HST Data 
Handbook ( 1 998 ), performed a sky background subtraction and 
used the synthetic zero points and transformation coefficients 
given by Holtzman et al. ( 1995b| ) to transform the measured 
fluxes to magnitudes in the Johnson U BVRI system. For the 
Ha narrow band filter, F656N, the continuum was subtracted 



as described in Section [2.1.l| ; the correction given by Eq. (|lj) 
was not necessary, as the Ha line lies outside the wavelength 
range covered by the F547M filter. The calibration of the Ha 
images was done following the prescriptions by Holtzman et 
al. (|1995b|). 



IRAF is distributed by the National Optical Astronomy Observa- 
tories, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 



2.2.2. Colour-magnitude diagrams 

The Hll region NGC2363 within Mkn 71 (Figure 1; right) 
was centered on the PCAVF1 chip (instrumental scale of 
07046/pixel), allowing the resolution of compact groups of 
stars with a mean linear separation of ~ 0.8 pc. We utilized the 
PCI data to derive CMDs in order to further constrain the re- 
sults of the spectral population synthesis analysis (Section 3.3) 
using the DAOPHOT II stellar photometry package under ESO 
MIDAS. A model for the point spread function (PSF), neces- 
sary for performing multiple-PSF fitting in crowded fields, was 
computed from several isolated bright stars in each frame. Be- 
cause of the severe crowding within some regions of NGC 2363 
as well as local background variations due to strong gaseous 
emission, we evaluated the photometric uncertainties and com- 
pleteness limits of the resulting point-source identification files, 
taking into account the error growth towards fainter luminosi- 
ties and the magnitude histograms (cf. Figure |5|). 

The measured fluxes were calibrated and transformed to the 
Johnson B and V filters as described in the previous section. 
The luminosity functions of the detected point sources (Figure 
||, left panel) suggest a reasonable completeness for sources 
brighter than « 26 mag in both B and V. Faintwards of the 
latter magnitude cutoff the photometric errors become domi- 
nant (see Figure Q). Close to the detection limit, an additional 
error source is introduced by the small-scale variations of the 
background due to gaseous emission as well as strong noise 
peaks which possibly are misinterpreted as faint stars (cf. the 
excess counts at low luminosities in the upper left panel of Fig- 
ure ||). We shall therefore consider the data reliable down to 
a limit of ~ 25 mag in B and V, close to the one Drissen et 
al. ( |1999| ) obtain for the same data set. This cutoff is sufficient 
for our purposes; from a simple comparison of data points and 
isochrones it becomes evident that, given the photometric er- 
rors and poor time resolution of the (B — V) CMD, no firm 
conclusions on the SF history can be drawn for ages il 50 Myr. 
A second CMD was derived in the same way for a region on 
the WF4 chip apparently free of any appreciable signatures of 
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Fig. 4. V vs. (B — V) colour-magnitude diagrams (CMDs) of Mkn71; overplotted are isochrones by Bertelli et al. (1994) for 
a metallicity of 1/20 Zq. Crosses to the left illustrate the mean errors in both axes as a function of the apparent V magnitude; 
the dashed lines indicate the estimated confidence limit for the B and V band, a CMD derived within the PCI chip of WFPC2, 
centered on the H II region NGC 2363; the isochrones correspond to ages of 4, 8, 12, 25, 50, 100 and 500 Myr. b CMD of the 
region of the main body of Mkn 7 1 indicated in Figure 1 (right panel, see label "CMD 2") where active star-forming regions are 
absent as derived from WF4 data. The isochrones correspond to ages of 25, 50 and 250 Myr. 



current SF (marked in Figure [Tj) to obtain information on the 
galaxy's underlying stellar population separately. From an in- 
spection of the right panels of Figure |^, and from the consid- 
erations referring to the CMD of NGC 2363 we estimate the 
confidence limit to «24 mag for both B and V. Either CMD 
is shown in Figure Q with typical photometric uncertainties for 
different magnitude intervals along with synthetic isoc hrone s 
for a metallicity of Z=0.001 adopted from Bertelli et al. (|l994|). 
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2.3. Spectroscopic observations and data reduction 

Spectrophotometric observations of Mkn 59 and Mkn 7 1 were 
obtained on March 17th and 18th 1994, using the Ritchey- 
Chretien spectrograph of the KPNO 4m Mayall telescope. We 
used a 2"x300" slit with the KPC-10A grating (316 lines 
mm -1 ) in first order, with a GG 385 order separation filter 
which cuts off all second order contamination for wavelengths 
blueward of 7400A. This instrumental setup gave a spatial scale 
along the slit of 0'.'69 pixel -1 , a scale perpendicular to the slit 
of 2.7 A pixel -1 , a spectral range of 3500-7500 A, and a spec- 
tral resolution of ~ 5 A. The seeing was l'/5 FWHM. Mkn 59 
was observed at two slit positions (cf. Figure |6j), both centered 
at the brightest star-forming region with PA. = 59° and PA. 
= 15° (close to the direction of the major axis). The total ex- 
posure time for the first orientation of the slit was 70 minutes 
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Fig. 5. Luminosity functions of the point sources included in 
the colour-magnitude diagrams (Figure |]) as obtained within 
the PCI- (left) and WF4 chip (right) of HST WFPC2. 
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Table 3. Emission Line Intensities in Mkn 59 



Ion region 1 region 2 







F(X)/F(H/3) 


J(A)//(H/3) 


F(X)/F(H/3) 




J(A)//(H/3) 


3727 


[On] 


1.0086±0.0018 


1.0879±0.0020 


1.9191±0.0832 


2 


.080 1±0 


.1049 


3750 


H12 


0.0204±0.0006 


0.0268±0.0009 










3770 


Hll 


0.0280±0.0006 


0.0349±0.0009 










3798 


H10 


0.0341±0.0006 


0.0415±0.0009 






... 




3835 


H9 


0.0498±0.0006 


0.0581 ±0.0009 










3868 


[Nem] 


0.4566±0.0011 


0.4872±0.0012 


0.3179±0.0503 





.3366±0 


.0584 


3889 


He I + H8 


0.1633±0.0008 


0.1786±0.0010 










3968 


[Nem] + H7 


0.2806±0.0009 


0.3014±0.0011 


0.1467±0.0309 





.2848±0 


.0698 


4026 


He I 


0.0110±0.0005 


0.0116±0.0005 










4069 


[Sn] 


0.0097±0.0006 


0.0102±0.0006 










4076 


[S II] 


0.0031 ±0.0036 


0.0033±0.0038 










4101 


H5 


0.2396±0.0008 


0.2552±0.0010 


0.1296±0.0246 





.2606±0 


.0598 


4340 


H 7 


0.4589±0.0011 


0.4766±0.0012 


0.3775±0.0316 





.4737±0 


.0461 


4363 


[Oiii] 


0.0890±0.0006 


0.0916±0.0006 


0.0560±0.0326 





.0550±0 


.0349 


4471 


He I 


0.0357±0.0005 


0.0365±0.0005 


0.0875±0.0320 





.0846±0 


.0339 


4658 


[III] 


0.0054±0.0004 


0.0054±0.0004 










4686 


Hell 


0.0131 ±0.0005 


0.0132±0.0005 










4713 


[ArlV] + He I 


0.0084±0.0004 


0.0085±0.0004 










4740 


[ArlV] 


0.0055±0.0004 


0.0056±0.0004 










4861 


H/3 


1.0000±0.0017 


1.0000±0.0017 


1.0000±0.0495 


1 


.0000±0 


.0554 


4921 


He I 


0.0082±0.0004 


0.0081 ±0.0004 










4959 


[OIII] 


2.0670±0.0030 


2.0503±0.0030 


1.4993±0.0672 


1 


.3588±0 


.0667 


5007 


[OIII] 


6.1585±0.0078 


6.0917±0.0078 


4.3242±0.1668 


3 


.8956±0 


.1645 


5199 


[Nl] 


0.0033±0.0004 


0.0033±0.0004 










5271 


[Fein] 


0.0019±0.0004 


0.0019±0.0003 










5518 


[CI III] 


0.0043±0.0004 


0.0041±0.0003 










5538 


[CI III] 


0.0026±0.0004 


0.0025±0.0004 










5876 


He I 


0.1138±0.0005 


0.1075±0.0005 


0.0854±0.0160 





.0698±0 


.0144 


6300 


[Ol] 


0.0210±0.0004 


0.0195±0.0004 










6312 


[S III] 


0.0209±0.0004 


0.0194±0.0003 










6364 


[Ol] 


0.0072±0.0003 


0.0066±0.0003 










6548 


[Nil] 


0.0242±0.0005 


0.0222±0.0005 










6563 


Ha 


3.0749±0.0041 


2.8214±0.0041 


3.6134±0.1384 


2 


.8062±0 


.1277 


6583 


[Nil] 


0.0540±0.0004 


0.0495±0.0004 


0.1800±0.0229 





.1377±0 


.0194 


6678 


He I 


0.0323±0.0004 


0.0295±0.0003 


0.0427±0.0137 





.0324±0 


.0114 


6717 


[Sn] 


0.1052±0.0005 


0.0959±0.0005 


0.2983±0.0266 





.2258±0 


.0225 


6731 


[S II] 


0.0785±0.0005 


0.0715±0.0004 


0.1597±0.0245 





.1207±0 


.0204 


7065 


He I 


0.0280±0.0003 


0.0252±0.0003 










7136 


[Arm] 


0.0862±0.0005 


0.0774±0.0004 


0.1331±0.0178 





.0974±0 


.0144 


7281 


He I 


0.0051±0.0003 


0.0046±0.0002 










7320 


[Oil] 


0.0186±0.0004 


0.0166±0.0004 










7330 


[On] 


0.0209±0.0005 


0.0187±0.0005 











C(Hf3) dex 
F(H/3) a 
EW(H(3) A 
EW(abs) A 



0.110±0.002 
15.08±0.02 
150.1±0.2 
0.35±0.04 



0.235±0.050 
0.28±0.01 
44±1 
3.9±0.5 



a in units of 10 ergs s cm 



and was broken up into four subexposures. At the second slit 
orientation a single 15 minutes exposure was taken. 

Mkn 71 was observed at a slit orientation of P. A. = 77° 
(see Figure []]), centered on the brightest star-forming region 
(NGC 2363 1). The total exposure time was 29 minutes, broken 



up into five subexposures. Three Kitt Peak IRS spectroscopic 
standard stars were observed during each night for flux calibra- 
tion. The airmasses during all observations were < 1.25, there- 
fore a correction for atmospheric dispersion was not necessary. 
Spectra of He-Ne-Ar comparison lamps were obtained before 
and after each observation to calibrate the wavelength scale. 
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Table 4. Emission Line Intensities in Mkn 59 



Ion 


reg 


ion 3 


F(\)/F(H/3) 


J(A)//(H/3) 


3727 [Oil] 


3.1730±0.1866 


2.9505±0.2018 


3868 [Nem] 


0.5371±0.0516 


0.4994±0.0528 


3968 [Ne in] + H7 


0.2177±0.0343 


0.2832±0.0603 


4101 HS 


0.1568±0.0426 


0.2377±0.0836 


4340 H7 


U.T^J JJ _1_VJ.VJ~J 1 


4RS3+0 0714 


4363 [Om] 


0.0594±0.0287 


0.0552±0.0287 


4861 H/3 


1.0000±0.0739 


1.0000±0.0848 


4959 [Om] 


1.3522±0.0939 


1.2574±0.0939 


5007 [Om] 


3.8327±0.2218 


3.5640±0.2219 


6300 [Oi] 






6563 Ha 


2.5651±0.1519 


2.4329±0.1693 


6583 [Nil] 


0.2177±0.0374 


0.2024±0.0378 


6717 [S 11] 


0.3142±0.0410 


0.2921±0.0418 


6731 [S 11] 


0.3037±0.0430 


0.2824±0.0437 


7136 [Arm] 






C(H/3) dex 


0.000±0.077 


F(H/3) a 


0.21±0.01 


EW(HP) A 


51.1±2.7 


EW(abs) A 


3.85±1.49 


a in units of 10~ 14 en 


>ss 1 cm 2 . 
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Fig. 6. B band image of Mkn 59; the two different slit positions 
and the regions investigated individually (see Section 2.3) are 
shown. The grey scale to the right displays the surface bright- 
ness levels in B mag arcsec -2 . 



The data reduction was done with the IRAF software pack- 
age. The two-dimensional spectra were bias-subtracted and 
flat-field corrected. Then the IDENTIFY, REIDENTIFY, FIT- 
COORD, TRANSFORM, BACKGROUND and CALIBRATE 
routines were used to perform the wavelength calibration, cor- 
rection for distortion and tilt, night sky subtraction and flux 
calibration for each frame. One-dimensional spectra were ex- 
tracted from the flux-calibrated two-dimensional spectra using 



the IRAF task APALL. For Mkn 59 we chose ten 7"x2" re- 
gions along the major axis of the galaxy (i.e the slit at RA. = 
15°), the locations of which are indicated in Figure || Six of 
these spectra are shown in Figure [7| From the slit at RA. = 59° 
the subspectra were extracted at each pixel row along the slit, 
so that the resulting areas were 0769 x 2" each. For Mkn 7 1 (see 
Figure [l]), the subspectra from regions 1 and 2 were extracted 
from 7"x2"regions, while for regions 3 and 4 areas of 14" x 2" 
and 28" x 2" were used to obtain a sufficiently high S/N . The 
extracted spectra are shown in Figure |[ 

The emission line intensities were measured utilizing a 
Gaussian profile fitting. All spectra were corrected for inter- 
stellar extinction, where the extinction coefficient C(H/3) was 
derived from the hydrogen B aim er decrement using the equa- 
tions given in Izotov et al. ( 1994 ) and the th eoretic al hydrogen 
emission line flux ratios from Brocklehurst ( 1971 1. The varia- 
tions of C(H/3) along different regions (Tables]f and \m arise 
mainly from observational uncertainties (imperfect focussing 
at the blue and red ends of the spectra), and from different com- 
binations of Balmer lines that could be reliably measured in 
each subspectrum to determine the local extinction. While each 
spectrum was corrected using its individual value of C(H/3), 
images were extinction-corrected adopting a uniform value of 
C(H/3)=0.09 and C( H/3)= 0.12 for Mkn 59 and Mkn 71, respec- 
tively (Guseva et al. |2000| ). 

The electron temperature X^Olll) was derived from the 
observed flux ratio [Om](A4959+A5007)/A4363. In cases 
where the [Om] A4363 emission line was not detected, we 
used the so-called upper branch of the Edmunds & Pagel 



(1984) calibration of the total oxygen emission line flux 
[Oll]A3727+[Om](A4959 + A5007) vs. the electron temper- 
ature to determine T e and derived the oxygen abundance fol- 



lowing van Zee et al. ( 1998a 1 



The observed (F{\)) and corrected (/(A)) emission line 
fluxes relative to the H/3 emission line fluxes for 3 regions in 
Mkn59 (slit RA. = 15°) are listed in Tables| andQ The tables 
contain only the regions where the emission line [Olll]A4363 
was detected at a S/N that allowed for a reliable flux mea- 
surement. Also listed are the extinction coefficient C(H/3), the 
observed flux of the H/3 emission line and its equivalent width 
EW(H{3) along with the equivalent width of hydrogen absorp- 
tion lines £W(abs). The line intensities of the brightest regions 
(1 and 2) of Mkn71 are presented in Izotov et al. (1997). 

Applying the electron temperature T e (Olll), ionic abun- 
dances of 2+ , Ne 2+ and Ar 3+ were derived. The tempera- 
ture T e (Oll) was inferred, according to Izotov et al. (1994, 
1996b, from the relation between T e (Oll) and T e (Olll) us- 



ing Hll region photoionization models by Stasihska (1990). 
From T e (0 II), the + and N + ionic abundances were deter- 
mined, while the intermediate value of the electron temperature 
T e (S III) served to derive the ionic abundances of Ar 2+ and S 2+ 



(Garnett |1992j). The [S ll]A6717/A6731 ratio was used to deter- 
mine the electron number density N e (S II). Total heavy element 
abundances were derived after co rrecting for unseen stages of 
i oniza tion following Izotov et al. ( 1994 1997 ) and Thuan et al. 
(|1995|). 
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The resulting ionic and heavy element abundances for the 
3 H II regions in Mkn 59 with measured [O III] A4363 are given 
in Table || along with the adopted ionization correction factors 
(ICF) while the heavy element abundances in the two brightest 



knots of Mkn 71 are given in Izotov et al. (1997) 



3. Results 

3.1. Structural properties of the host galaxies 

Recent studies have established that the underlying host galaxy 
of iE/nE BCDs does systematically differ with respect to its 
central surface brightness [ie o and exponential scale length a 



from other classes of dwa rf galaxies suc h as d ls and dEs ( P96b , 
Patterson & Thuan 1996 Marlowe et al 
Marlowe et al 



Salzer & Norton 1998 



1997 



Papaderos 



1998, 



1999). This structural 



dichotomy is evident from Figure ^| (adopted from Papaderos 



1998) showing that the central surface brightness and exponen- 



tial scale length of the LSB component of a BCD with an ab- 
solute B magnitude Me ^-16 mag are respectively ^1.5 mag 
brighter and a factor of ~ 2 smaller than in a typical dl/dE 
of equal luminosity. Although there is no sharp limit, a gap 
around He,q ~22 Bmag/d" separates the host galaxies of 
iE/nE BCDs from other classes of dwarf galaxies. The same 
diagram shows that more luminous BCDs (i.e. systems with 
a host galaxy brighter than -16 mag in the B band) follow 
the same trend populating systematically different areas in the 
He,o~Me and log(a£)-ME-planes than dls and dEs. In Fig- 
ure ^ (left panel) we show with open squares the positions of 
the LSB components of the LLC BCDs studied here (values 
derived in the B band) alon g with the ones of the iI,C BCDs 
SBS 1415+437 (Thuan et al. |l999|), To l 1214-277 (Fricke et al. 
2000| ) and UM417 (Cairos et al. ^000| ). In the right panel of the 
same figure we include two further iI,C BCDs for which mea- 
surements of the exponential scale length of their LSB com- 
ponents are availa ble, U M 133 (Telles & Terlevich 1997 ) and 
Mkn 1328 (James |1994| ). It may be seen that the iI,C BCDs fit 
into the gap between typical iE/nE BCDs and dls/dEs except 
for SBS 1415+437 and Tol 1214-277 which are indistinguish- 
able from typical BCDs in the he,q-Me and \og(otE)-ME pa- 
rameter space. 

The uncertain distance to Mkn 59 (cf. Section [l]) does not 
affect the results stated above. As can be seen in both panels 
of Figure |9j the vector illustrating a shift to the data points due 
to a change of the Hubble constant is approximately parallel to 
the sequences of data points for iE/nE BCDs and other types 
of dwarf galaxies. For a different distance, Mkn 59 therefore 
remains within the gap between compact BCDs and dls/dEs. 

3.2. Heavy element and line intensity distributions 



In Figure 10 we show the spatial intensity distributions of the 



Ha and H/3 lines along with that of the continuum adjacent to 
H/3, as determined from the one-dimensional spectra extracted 
along the slit at PA. = 59° of Mkn 59 (cf. Figure 6). It is 
obvious that all three distributions are spatially coincident. 



This is also the case for other nebular emission lines, including 
He I, [Nil], [Oil], [Om] and [Nelll], which show the same 
spatial distribution as the Ha and H/3 emission lines. 

Wolf-Rayet stars were detected in region No. 1 of both, 
Mkn 59 and Mkn 71. For Mkn 71, 5 late nitrogen WR stars 
(WN stars) and 2 early carbon WR stars (WC stars) stars were 
found. The number of WC stars in Mkn 59 was determined to 



be 13 and that of WN stars to be 40 (Guseva et al. 200C , values 



adapted to the distances used here). The spatial flux distribu- 
tions of the blue A4650 and red A5808 WR bumps are shown 
in the right panel of Figure [l0| along with the one of the nebu- 
lar He II A4686 emission line. The red bump results primarily 
from WC stars, the blue bump from WN stars. We remark that 
the maxima of the H/3 and WR emission are shifted by 1 pixel 
{Of! 69) along the slit at PA. = 59° suggesting that the locations 
of O and WR stars may slightly differ. 

In all one-dimensional spectra of Mkn 59 at PA. = 59° the 
[O III] A4363 emission line is detected which allows a reliable 
determination of the electron temperatures and element abun- 
dances. Despite the large spatial variations of the oxygen line 
intensities and their ratios, the oxygen abundance is practically 
constant over the starburst region within the errors (Figure pTj ). 
The abundances of other heavy elements (neon, nitrogen, ar- 
gon, sulphur, iron and chlorine) show the same, nearly constant 
spatial distribution along the slit with PA. = 59° through the 
center of the giant H II region. In spite of the small abundance 
variations within the central part, in most cases < 0.2 dex, we 
remark that a weak gradient is present for the abundances of 
all heavy elements along the slit from -2" to +3", i.e. on scales 
of ~ 260 pc. This suggests a local heavy element enrichment 
and may be related to the possible displacement of WR stars 
relative to O stars (see above). 

In Figure [j~2| we show the spatial distribution of the oxy- 
gen abundance in Mkn 59 along the major axis with PA. = 
15°. Three regions where the electron temperature was derived 
from the [O III] A4363/(A4959 + A5007) flux ratio are indi- 
cated by stars, while open circles mark the regions where the 
[O III] A4363 line was not detected, thus oxygen abundances 
were derived following van Zee et al. (1998a). The latter abun- 
dances appear systematically higher compared to those where 
the electron temperature could be directly constrained utiliz- 
ing the [O III] A4363 line. However, the differences are still 
comparable to the intrinsic uncertain ty of th e order of ±0.2 in 
log(OZH) ascribed by van Zee et al. ( 1998a ) to their empirical 
calibration. Taking this fact into account, as well as the weak- 
ness of the [O III] A4363 line in the outer parts of the galaxy, 
we shall consider the oxygen abundance as constant over the 
main body of Mkn 59 with an average value of 12+log(0/H) = 
8.0. 

The same analysis for Mkn 71 yields for the two bright- 
est regions (regions 1 and 2; Figure [j], right), in which the 
[Olll]A4363 line is observed, heavy element abundances of 
12+log(0/H) = 7.79 and 7.77, respectively. The oxygen abun- 
dance in the regions 3 and 4 derived from the empirical cali- 
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Table 5. Heavy Element Abundances in Mkn 59 



Property 


region 1 


region 2 


region 3 


T e (OlIl)(K) 


13,469±38 


13,069±3,299 


13,558±2,910 


T e (Oll)(K) 


12,979±34 


12,762±3,008 


13,027±2,626 


T e (SlIl)(K) 


13,224±32 


12,916±2,738 


13,293±2,415 


N e (S Il)(cm- 3 ) 


79±10 


10±10 


536±561 


O+/H+(xl0 5 ) 


1.502±0.012 


3.001±2.164 


4.239±2.571 


O++/H+(xl0 5 ) 


8.757±0.069 


6.152±4.365 


5.095±3.013 


O +3 /H+(xl0 5 ) 








O/H(xl0 5 ) 


10.26±0.070 


9.154±4.872 


9.334±3.961 


12 + log(0/H) 


8.011±0.003 


7.962±0.231 


7.970±0.184 


[0/H] a 


-0.919±0.003 


-0.968±0.231 


-0.960±0.184 


N+/H+(xl0 7 ) 


4.902±0.040 


14.26± 10.28 


19.98±12.12 


ICF(N) 


6.83 


3.05 


2.20 


log(N/0) 


-1.486±0.006 


-1.323±0.550 


-1.327±0.455 


Ne ++ /H+(xl0 5 ) 


1.639±0.014 


1.243±0.931 


1.646±1.016 


ICF(Ne) 


1.17 


1.49 


1.83 


log(Ne/0) 


-0.728±0.006 


-0.695±0.461 


-0.491 ±0.374 


S+/H+(xl0 7 ) 


2.164±0.013 


4.634±1.900 


7.776±2.782 


S ++ /H + (xl0 7 ) 


14.85±0.281 






ICF(S) 


1.89 






log(S/0) 


-1.504±0.006 






Ar ++ /H+(xl0 7 ) 


3.497±0.023 


4.652±1.738 




Ar +3 /H+(xl0 7 ) 


1.049±0.080 






ICF(Ar) 


1.02 


1.53 




log(Ar/0) 


-2.346±0.008 


-2.108±0.266 




Fe++/H+(xl0 7 ) 


1.505±0.123 






ICF(Fe) 


8.54 






log(Fe/0) 


-1.902±0.012 







a [X] = logX - logXo . Solar abundances are from Anders & Grevesse ( 1989| ) 



bration method described in van Zee et al. (1998a) is larger by 
-0.3 dex. 



3.3. Population synthesis models 

For an analysis of the stellar content of the young ionizing 
clusters, we used the spectral e nergy distributions (SEDs) cal- 
culated by Schaerer & Vacca ( 1998 ) for heavy element mass 
fractions and ages in the range between Z = 0.001-0.02 and t 
= 0.1 - 10 Myr, respectively. For ages t > 10 Myr, we calcu- 
lated a grid of SEDs for stellar populations with ages ranging 
from 10 Myr to 10 Gyr in time steps of Alogt[yr] =0.1 and 
a heavy element mass fraction of Z = 0.002, using isochrones 
from Bertelli et al. (1994) and the compilation of stellar at- 
mosphere models from Lejeune et al. (1998). An initial mass 
function (IMF) with a Salpeter slope (2.35) and lower and up- 
per mass limits of 0.6 Mq and 120 M Q was adopted. 

To study the age of the stellar populations in the galax- 
ies and to compare the results of the spectral analysis with 



broad band photometric data, it is necessary to take into ac- 
count both the stellar and ionized gaseous emission in the spec- 
tra. For this purpose, the stellar SED has to be separated from 
the gaseous emission following the procedure described by Gu- 
seva et al. ( |1998| ) and Papaderos et al. ( |1998| ). We added stellar 
SEDs calculated for instantaneous bursts with different ages 
("single stellar populations"; SSPs) to the observed gaseous 
emission SED to match the total SED. The contribution of the 
gaseous emission was scaled to the stellar emission by the ra- 
tio of the observed equivalent width of the H/3 emission line 
to the equivalent width of H/3 expected for pure gaseous emis- 
sion. To calculate the gaseous continuum SED at each region 
along the slit, the observed H/3 flux and the electron temper- 
atures were derived from the respective spectrum. The contri- 
butions of bound-free, free-free, and two-photon processes to 
the continuum emission were then c alcul ated fo r the spectral 
range from to 5 /im (Ferland |1980| , Aller |1984| ). 
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Fig. 7. Spectra of different regions of Mkn 59 extracted along the major axis slit at P.A. = 15° (thin lines). Overplotted are 
spectral energy distributions (SEDs) of model composite stellar populations and gaseous SEDs. The SEDs of regions No. 3, 7 
and 10 are redder than those of the other regions and refer to the underlying stellar component. It is obvious that the contribution 
of ionized gas to the observed SED is much larger in region 1 than in other regions further to the northeastern direction of the 
main body of the BCD. 



3.4. The ages of the stellar populations 

The number of O stars within the giant H II region complexes at 
the southeastern tips of Mkn 59 and Mkn 7 1 was derived from 
the H/3 flux to w 4740 and 2 010, r espectively, following the 
prescriptions by Guseva et al. ( 200C ). As stated in Section l_2, 



WR stars were detected in both BCDs, which, together with the 
large number of O stars, contribute strongly to the total SED. 
Taking into account the average metallicities, the ages of the 
brightest H II regions in which WR stars were detected cannot 
exceed 4 to 5 Myr (e.g. Schaerer & Vacca 1998). 
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Fig.8. Spectra of the regions in Mkn71, extracted from the slit at P.A. = 77° (thin lines). Overplotted are spectral energy 
distributions (SEDs) of model composite stellar populations and gaseous SEDs. The SED of the region No. 4 is redder than those 
of other regions and matches the underlying stellar component. Note the large contribution of the ionized gas in region 1 . 



In Figs. ^ and |8| the modelled stellar SEDs, the gaseous 
SEDs and their coadded spectral distributions (marked as "to- 
tal") are shown superposed on the observed spectra for different 
regions of the BCDs. For nearly all spectra, a good agreement 
is achieved between the observed and modelled SEDs. Only in 
the case of region 2 of Mkn 71 (see Figure [j], right and Figure 
|8|), the agreement is not as good due probably to an imperfectly 
focussed telescope. 

Assuming that the observed gas emission is not caused by 
local SF, but due to photoionization or shock ionization by SF 
outside the observed region, the observed SEDs are in some 
cases reproducible with one SSP of intermediate age (few 10 8 
yr) only; mostly, however, reproducing the observed SEDs re- 
quired a superposition of SSPs with different ages and mass 
fractions. Tables ^ and ^ list the ages and relative mass frac- 
tions of the old and young SSPs which provide the best fits 
to the observed SEDs at different slit positions. The models 
imply respective ages of 4 and 2.9 Myr for the stellar popula- 
tion formed in the bright H II region complexes in Mkn 59 and 
Mkn 71 (slit regions 1) while for other slit positions the age of 
the younger stellar continuum was inferred to few 10 Myr. The 
ages of the old SSPs could be constrained to ^ 2 Gyr. Within 



the uncertainties of the methods applied, an upper limit of ~ 
3-4 Gyr is possible (see the following subsections), whereas 
even higher ages appear untenable. One has, however, to keep 
in mind that present spectrophotometric dating methods cannot 
definitely rule out the presence of a small fraction of even older 
stars, which due to their high M/L would barely contribute to 
the SED. 

3.4.1. Uncertainties of the population synthesis models 

The goodness of the spectral fits was found to sensitively de- 
pend on the mass fraction and ages of the adopted young and 
old SSPs. For instance, for region 1 of Mkn 59, increasing the 
age of the young SSP from 4 to 5 Myr yielded a satisfactory 
fit, but a 6 Myr old SSP resulted in systematic residuals be- 
tween the observed and modelled SED. The age of the old SSPs 
could be constrained with an accuracy of w 1 Gyr at slit regions 
where young stellar sources provide a minor contribution to the 
light (regions 3, 7 and 10 in Mkn59 and region 3 in Mkn71). 
In e.g. region 7 of Mkn 59, old SSPs up to 2 Gyr could still 
reproduce the observed SED, while already a 3 Gyr old SSP 
produced obvious residuals to the observed spectrum. A com- 
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Fig. 9. B band structural properties of the underly ing ex ponential components for different types of dwarf and low surface 
brightness (LSB) galaxies, compiled in Papad eros (1998| ). The positions of the com etary BCDs Mkn59, Mkn71, Mknl328, 
UM 133, UM417, SBS 1415+437 (Thuan et al. |1999| ) and I Tol 1214-277 (Fricke et al. |2000| ; cf. Section this work) are indi- 
cated by big squares, left Extrapolated central B surface brightness (iE,o vs. absolute B magnitude Me of the LS B components, 
right E xpone ntial scale length a vs. Me- Boxes with arrows to the right mark UM 133 (Telles & Terlevich 1997 1 and Mkn 1328 
(James 1994) and give upper limits for Me- Data points for iE/nE BCDs and other classes of dwarf galaxies are taken from 
Drinkwater & Hardy (|1991[ ), ^964 Papaderos ( |1998|) Marlowe et al. ( |1997|) , Noeske ( |1999| ), Noes ke et a l. ( |l999a| > and B inggel i 



& Cameron (1991, 1993), Caldwell & Bothun (1987), Bothun et al. (1991), Patterson & Thuan (1996), Vigroux et al. (1986) 



Hopp & Schulte-Ladbeck (1991) and Carignan & Beaulieu (1989), respectively. 



parison of the observed and modelled EW(Hj3) yielded further 
constraints to resolve ambiguities between different SED solu- 
tions. The assumption of high ages for the old SSP requires 
an increased contribution by the young SSP to reproduce the 
blue continuum, which in turn results in too high EW(H(3) as 
compared to the observed value. 

Given that the metallicities of the stellar populations in 
BCD s may be low er tha n those of the ioniz ed ga s (Calzetti 
1997| , Guseva et al. |l998j Mas-Hesse & Kunth |l999| ), we com- 



puted a set of models varying the metallicity of the SSPs around 
Z =0.002 by 0.5 dex, the smallest stepsize our model libraries 
allowed. Models computed on the assumption of a metallicity 
other than Z =0.002 failed to adequately reproduce the ob- 
served spectrum, suggesting that the metallicity of the ionized 
gas is similar to that of the stellar population. 

Alternatively, to constrain the age of the underlying galax- 
ies, not only an instantaneous burst was considered, but also 
an extended episode of moderate SF at a constant rate. This 
way, the presumably complex formation histories of the galax- 
ies were bracketed between two limiting cases. 



In the same way as described in Section 3.3, a grid of model 



SEDs was calculated for stellar populations that were formed 
in extended SF episodes which started and ended at different 
times. The best-fitting solutions for e.g. Mkn 71 correspond to 
a formation episode of the underlying galaxy which started ~ 



2 Gyr ago and ended ~ 10 Myr ago. Within the uncertainties, a 
SF episode that began up to 3 Gyr ago is possible; higher ages 
of the underlying galaxies would require a very recent (< 10 
Myr ago) termination of their formation episode to reproduce 
the observed SEDs. Such a recently ongoing SF activity in the 
underlying galaxy in Mkn 71 can be rejected from the colour- 
magnitude diagrams (Figure Q); in the left panel (a), obtained 
at the starburst region NGC 2363, stars are absent between the 
actual burst (younger than ^10 Myr) and the last SF episode of 
the underlying galaxy which seems to have occurred more than 
~ 20 Myr ago (note the red supergiants at the tip of the 25 Myr 
isochrone). Also in the right panel (b), no stars much younger 
than ~ 20 Myr are present. This indicates that, in the galaxy 
underlying the starburst, no significant SF occured during the 
last > 20 Myr, in agreement with the ground-based CMDs by 
Aparicio et al. (1995). A similar analysis for Mkn 59 indicates, 
within the unvertainties, an upper age limit of ~ 4 Gyr. 

3.4.2. Consistency with imaging data; further implications 

The EW(U/3) of region 1 in Mkn 59 and Mkn 71 amount to 
150A (Table g) and 316A(Izotov et al. 1997), respectively. For 
these regions the broad band colours calculated from the syn- 
thetic SEDs, including and ignoring gaseous emission (Tables 
H and [7]) differ from each other by up to 0.5 mag. This demon- 
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Fig. 10. left Spatial distribution of the continuum flux near H/3 (dashed line) and of the fluxes of the H/3 and Ha emission lines 
(thin and thick solid lines, respectively) in Mkn 59 in the direction of P. A. = 59°. The fluxes are normalized to the values at the 
brightest pixel, right Spatial distribution of the WR broad bump fluxes at A4650 (WN bump, thick line), A5808 (WC bump, 
thin line) and the nebular He II A4686 emission line flux (dashed line) in Mkn 59 along the P.A. = 59° direction. Note that the 
maximum of the WN bump appears to be slightly shifted relative to the maximum of the H/3 emission line, continuum and 
nebular flux distribution. 
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Fig. 11. Spatial distributions of the heavy element abundances in Mkn 59 along the slit at P.A. = 59° 



strates that a reliable dating of stellar populations within a star- 
burst environment requires a correction for the colour shift in- 
duced by ample gaseous emission. 

The consistency of the observed broad band colours with 
the spectral population synthesis results was checked by ex- 
tracting the areas covered by the spectrograph slits from the 
colour maps of Mkn 59 and Mkn 71 (cf. Section |2.1.1[ ). Fig- 
ure [l3| shows the modelled and observed colours along the slits 
after the latter were transformed to the Johnson-Cousins sys- 
tem (Bessell 1990| ). Generally, the observed and synthesized 
colours agree well with each other (within la uncertainties), 
given that i) the synthesized colours are averaged values over 
regions covering several arcsec length along the slit, and ii) a 
minor displacement between the positions of a region in the 
colour map and of the slit may significantly alter the observed 
colours (see Figure || and the 1 cr-variations of the observed 



colours in Figure [jj|), as the H II regions are typically very 
compact. 

Transforming the colours of the LSB component as derived 
from the colour profiles (Figure [|) to the Johnson-Cousins sys- 
tem yields a B — R <>0.6 for Mkn 59 and Mkn71. This is in 
good agreement with the upper limits of the colours derived 
from the spectra (cf. Tables ^ and fy. The constancy of the 
colour profiles for photometric radii where the starburst light 
contribution becomes negligible (i.e. for R* <; P25) indicates 
that in the underlying galaxies population gradients are absent. 
Therefore, the upper age limits as being derived for either BCD 
along the spectrograph slits can be considered to be valid for 
the entire underlying galaxy. 

4. Discussion 
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Table 6. Synthetic colours of a mixture of young and old populations in Mkn 59 



region 


log* 


relative mass 


logt 


relative mass 


(B-V) 


(B-R) 


(R-I) 


(V-K) 








a) Stellar and gaseous emission 








1 


6.6 


0.7 


9.0 


0.3 


0.17 


-0.01 


-0.32 


0.12 


2 


6.9 


0.6 


9.3 


0.4 


0.03 


0.09 


-0.03 


0.64 


3 


8.7 


1.0 






0.29 


0.52 


0.28 


1.93 


4 


7.4 


0.1 


9.3 


0.9 


0.21 


0.42 


0.25 


1.59 


5 


7.4 


0.1 


9.3 


0.8 


0.28 


0.44 


0.16 


1.46 


6 


7.0 


0.3 


9.3 


0.7 


0.04 


0.21 


0.24 


1.39 


7 


8.7 


1.0 






0.28 


0.58 


0.30 


2.03 


8 


7.4 


1.0 






0.05 


0.24 


0.27 


1.54 


9 


7.0 


0.1 


9.3 


0.9 


0.11 


0.30 


0.28 


1.51 


10 


8.3 


1.0 






0.15 


0.33 


0.32 


1.90 










b) Stellar emission 








1 


6.6 


0.7 


9.0 


0.3 


-0.15 


-0.17 


0.00 


0.20 


2 


6.9 


0.6 


9.3 


0.4 


-0.04 


0.02 


0.11 


0.70 


3 


8.7 


1.0 






0.28 


0.54 


0.39 


2.08 


4 


7.4 


0.1 


9.3 


0.9 


0.18 


0.43 


0.33 


1.70 


5 


7.2 


0.2 


9.3 


0.8 


-0.01 


0.11 


0.19 


1.22 


7 


8.7 


1.0 






0.28 


0.54 


0.39 


2.08 



Table 7. Synthetic colours of a mixture of young and old populations in Mkn 7 1 



region 


logt 


relative mass 


logt 


relative mass (73 — V) 


(B-R) 


(R-I) 


(V-K) 








a) Stellar and gaseous emission 








1 


6.46 


1.0 




0.26 


0.04 


-0.58 


0.07 


2 


6.6 


0.15 


9.3 


0.85 0.03 


0.13 


-0.14 


0.62 


3 


6.6 


0.02 


9.0 


0.98 0.30 


0.57 


0.40 


2.09 


4 


8.7 


1.0 




0.28 

b) Stellar emission 


0.54 


0.39 


2.08 


1 


6.4 


1.0 




-0.19 


-0.25 


-0.08 


-0.42 


2 


6.6 


0.15 


9.3 


0.85 -0.07 


-0.03 


0.09 


0.67 
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Fig. 12. Spatial distribution of the heavy element abundances 
in Mkn 59 along the slit at P. A. = 15°. Stars mark the regions 
where the [O III] A4363 emission line was detected, open cir- 
cles correspond to regions where the oxygen abund ance w as 
derived from the empirical relation by van Zee et al. (1998a). 



4.1. The metallicity distribution 

Compared to BCDs in general, which show a metallicity dis- 
tribution that peaks around w 1/10 Zq and steeply decays to- 



wards lower metallicities (Kunth & Sargent |1986| ), the metal 
abundances of Mkn59 (1/8%)) and Mkn71 (1/14Z Q ) are 
not exceptional. The nearly constant metallicity along the ma- 
jor and minor axes of Mkn 59, as well as along the minor axis 
of Mk n 7 1 (Figs, [fl] and [l|; see also Gonzalez-Delgado et al. 



1994), suggests that large-scale mixing processes in the ISM 
of both BCDs were at work on time scales of a few 10 6 yr. 
Otherwise, one would expect measurable metallicity enhance- 
ments in the vicinity of regions of ongoing or recent SF. Fur- 
thermore, the transport of metal-enriched warm and hot gas on 
scales of up to ~ 1 kpc from the starburst region, an observa- 
tional signature of which is the formation of supergiant shells, 
is commonly seen in th e ISM of star-forming dwa rf gala xies 



(cf. e.g. M arlowe et al 
mans et al 



1995| , Hun ter & Gallagher |1997 



Bo- 



1997| , Brinks & Walte r |1998| Papaderos & Fricke 
1998| , Strickland & Stevens [l999l cf. also Figure |, this work). 



The dilution of heavy elements in the vicinity of starburst re- 
gions may even be powered by galactic outflows. These were 
predicted b y num erical models to develop within dwarf galax- 
ies (Vade r fl986t De Young & Galla gher |1990[ De Young & 



1994, Mac Low & Ferrara 



1999) and proved neces- 



Heckman 

sary to account for the much lower metallicities observed in 
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Fig. 13. Observed colour distributions along the slits extracted 
from colour maps, compared to the colours computed from 
population synthesis solutions (black dots) for the analyzed 
subspectra. All colours are derived on the assumption of a uni- 
form extinction and are transformed to the standard Johnson- 



Cousi ns UB VRI system specified by Bessell ( 199C ), Bessell 
et al. ( |1998| ). top: Mkn59, P.A. = 15°, bottom: Mkn71, P.A. = 
77°. Vertical error bars of the synthetic colours include cumula- 
tive calibration uncertainties of the colour maps and population 
synthesis models, and of the transformation between the filter 
systems; horizontal error bars give the length of regions along 
the slit from which the individual subspectra were extracted. 
The alignment uncertainty between the observed and synthetic 
colors in direction of the slit (i.e. along the abscissa) is ps 5" . 
The error bars of the observed colour distributions give the la 
colour change for a 2" shift of the extracted colour map region 
perpendicular to the slit direction. 



BCDs with respect to those predicted from closed-box evo- 
lutionary synthesis models (Kriiger 1992 Lisenfeld & Ferrara 



1998). The development of such a large-scale perturbation of 
the ambient gaseous component does not appear unreasonable 
given the high-velocity (~ 10 3 km sec -1 ) gaseo us motions re- 
vealed spectroscopically by Izotov et al. ( 1996 ) and Roy et al. 
( 1991 see also Figure |[ this work) in Mkn59 and Mkn71, 
respectively. 

The small variations of the heavy element abund ance s 
along the giant Hll region complex of Mkn59 (Section 3.2), 
as well as the possible spatial offset of the O- and WR stars, 
may be explained by a relocation of SF processes within a time 
span of few Myr. Propagating SF on a linear scale of ~ 400 pc 
wit hin th e last 10 Myr has also been suggested by Drissen et 
al. (1999), to account for the age differences of the young star 
clusters NGC 2363 I and II within Mkn 7 1 . 



4.2. Morphology vs. structural properties and age 

The two iI,C BCDs investigated here appear similar with re- 
spect to the age and structure of their exponentially distributed 



host galaxies. This raises the question of whether the entire 
class of iI,C BCDs shares, besides a morphological resem- 
blance, a set of common physical properties, such as an inter- 
mediate age and structural properties bridging the gap between 
iE/nE BCDs and dls/dEs. A literature search for data on other 
nearby cometary BCDs does indeed provide some support to 
this hypothesis. 

Mkn 1328 /VCC 1374: NIR surface photometry (James 
1994) yields an exponential scale length of a=0.88 kpc for its 
LSB component. With its integrated absolute B-magnitude of 
-16 mag, and taking into account that a burst raises the B lu- 
minosity of a BCD by typically ~ 0.75 mag (P96b, Salzer & 



Norton 1998), Mkn 1328 resides presumably in between the 
parameter spaces populated by BCDs and dls in the log(as)- 
Mg-plane (cf. Figure 9). From the integral (B — H) colour of 
« 0.6 mag, the age of its host galaxy may be estimated to 
6 Gyr, following the predictions by Kriiger ( 1992j ) and Kriiger 
et al. (1995) when a burst parameter b < 0.01 is adopted. The 



metal licity was determined to 1/8 Zq (Kinman & Davidson 
198TJ). 

UM133: For the observed metallicity (1/17 Zq, Telles 



1995 ), the integral (R — I) colour of 0.4 mag derived by Telles 
& Terlevich (1997) is consistent with an age of ~ 4 Gyr. With 
My=-18.25, and an exponential scale length of a=1.39 kpc for 



its LSB component (Telles & Terlevich 1997), this object falls 
also into the gap between BCDs and dls. 

UM417 (Zk. 1/13 Z Q ; Campos-Aguilar at al. |1993| ). Sur- 
face photometry of the LSB component yields an intermedi- 
ately compact structure, too (ag = 0.45 kpc, Me = -14.4 B 
mag and p E .o = 22.4 B mag/d"; Cairos et al. 2000| ). The same 
authors deduce a (B — V) color for the LSB component of Si 
0.5 mag, compatible with an age of few Gyr. 



SBS 1415+437 (Thuan, Izotov & Foltz |1999|: With a prob- 
able age of <> 100 Myr and an extremely low metallicity 
(1/21 Zq), this galaxy is a young galaxy candidate. Its LSB 
component, however, does not fit within the structural gap be- 
tween BCDs and dls; with fj, E .o=2l.O Vmag/D", OE=0.30kpc 
and M E =-14.95 V mag, it is rather comparable to a compact 
BCD. On the other hand, as this object apparently still under- 
goes its first major episode of SF, it is likely that its stellar LSB 
component has not yet been fully built but, as suggested by 
the age gradient along its major axis, delineates the trail along 
which SF has occurred. 

Toll214-277: Similar to SBS 1415+437, this is an ex- 
treme ly low -metallicity iI,C BCD (Z » Zq /23, Terlevich 
et al. 1991). Recent surface photometry studies with the VLT 
(Fricke et al. |2000[ ) yield a E w 0.48 kpc and M E w -16 B 
mag. The average colours of the LSB component, (U — B) <, 
-0.4 mag and (B — R) w +0.3 mag, suggest an unevolved 
galaxy with an age of few 10 8 yr. 

The examples given above call for a further investigation of 
the following two hypotheses: 

(i) Cometary BCDs are relatively young objects, with ages not 
exceeding a few Gyr, thus systematically younger than the ma- 
jority of BCDs classified iE/nE. 
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(ii) Except for the candidate young extremely metal-deficient 
galaxies (r ~ 10 s yr) among them, cometary BCDs show 
structural properties of their host galaxies being typically in- 
termediate between those of iE/nE BCDs and dls/dEs. 

Provided that the above hypotheses will be strengthened by 
an investigation of a larger sample of nearby iI,C BCDs, one 
may expect the iI,C morphology to occur more frequently in 
dwarf objects at higher redshifts. Indeed, examples of high-z 
galaxies displaying a comet-like morpho logy have frequently 
been reported (cf. e.g. Dickinson 1996 for a galaxy cluster 
at 0=1.15). Furthermore, a few medium-redshift objects with 
comet ary m orphology were included in the sample of Guzman 
et al. ( |1998| , Figure 1), notably the Faint Blue Galaxy (FBG) 
HERC 13088, at a redshift z=0.436. With a total luminosity 
of Afg=-21.5 mag and an exponential scale length a w 2.8 
kpc, this galaxy appears an upscaled version of the nearby 
iI,C BCDs studied here. It is worth noting that the rest-frame 
(B - V) colours of the LSB host of the latter FBG, ~ 0.35 
mag, also suggest an age of ^ 5 Gyr. 

These findings provide support to the hypotheses proposed 
in Section 1, namely that a strong off-center burst ontop a 
dwarf LSB component as observed in iI,C BCDs is primar- 
ily not a stochastic event but occurs in systems considerably 
younger than iE/nE BCDs, probably differing from the latter 
with regard to their structural properties. These results empha- 
size the need for a detailed investigation of the processes lead- 
ing to the development of iI,C morphology in dwarf galaxies. 
We shall briefly remark on that issue in the next section. 

4.3. An elongated structure of high gas density? 

In the majority of BCDs, the centrally concentrated H II com- 
plexes where massive SF takes place are typically found to al- 
most coincide with maxima in the surfa ce density of an ex- 
tended Hi envelope (Taylor et al. 1994, Simpson 1995, van 



with a probably primordial chemical composition (Thuan & 



Zee et al. |1998bJ). Since a high gas density is necessary to sus- 
tain SF, such a condition is expected to be fulfilled near the 
end of the stellar body of an iI,C BCD. Age gradients of the 
star-forming regions along the stellar body's major axis seem 
to be frequent attrib utes of such a syst em (Barth et al. 1994 ; 
Aparicio et al. 1995 ; Thuan et al. 1999 ). These may be plausi- 
bly explained as due to propagation of SF activity (cf. Thuan et 
al. 1987) along the major axis of an elongated H I body. 

Given these observational results, it seems that the H I halos 
of iI,C BCDs have elongated central concentrations with major 
axes that coincide with those of the stellar bodies. Indeed, low- 



resolution H I maps of Mkn 59 ( Wilcots et al. 1996 ) and Mkn 7 1 



(Wevers et al. 1986) show the stellar components to follow pro- 



late central concentrations within elongated H I clouds. A sim- 
ilar distribution of the gaseous halo does not seem rare among 



magellanic irregulars (Wilcots et al. 1996) and has been inter 



ferometrically mapped also in BCDs in different evolutionary 



stages (see e.g. II Zw40, van Zee et al. 1998b, andlZw 18, van 
Zee et al. |l998c|). Moreover, the young BCD SBS 0335-052 
(Izotov et al. 1990) is forming within an elongated gas cloud 



Izotov 1997, Lipovetsky et al. 1999) 



It might be argued that stellar bodies and H I concentrations 
of iI,C BCDs merely represent edge-on disks, as suggested by 
the slow rotation found for the H I envelopes of Mkn 59 and 
Mkn 71 as well as for the ionized gas of SBS 1415+437 (Thuan 
et al. 1999| ). The generally exponential surface brightness pro- 
files of their stellar LSB hosts (cf. Fi gure p| ) do not serve as 
a proof of a disk structure (Freeman 1970), since exponen- 



tial SBPs are also observed for e.g. small spheroidal systems 



(Binggeli et al. 1984) and bars in late-type spirals (Elmegreen 
et al. 1996| ). However, the assumption of edge-on disks would 
demand the presence of numerous less-inclined iI,C BCDs, i.e. 
relatively blue disks with a single starburst in their outskirts, in 
contradiction to morphological studies of dwarf galaxies. 

4.3.1. iI,C systems vs. old iE/nE BCDs 

Considering that only little data is available for BCDs with 
morphological types other than iE/nE, the construction of a 
tentative age-morphology sequence for gas-rich dwarf galax- 
ies appears premature. The evolution of the different types of 
dwarf galaxies as well as their possible relations among each 
other are still sketchy, especially the triggering and effects of 
starbursts as well as the role of the Dark Matter (see Thuan 
Davies & Phil lipps |1988|, f96b| , Meurer |1998|, S waters 



1985 



1998 



Pazetal. 1999) 



Marlowe et al. 1997, 1999; van Zee et al. 1998b, Gil de 



However, we argue that their relatively low ages would 
make iI,C BCDs a possible link between the extremely young 
(t ~ 100 Myr) dwarf galaxy candidates, most of which be- 
long to the iO class, and the evolved iE/nE BCDs, for which a 
lower age limit of 2 Gyr has been invariably derived (Kriiger & 
Fritze-v.Alvensleben 1994 Noeske |l99^ ). In the young BCDs 
investigated thus far, a propagation of the SF activity has been 
observed (Papaderos et al. 1998, Izotov et al. 1999| , Thuan et 



al. 



1999). Thereby, it is conceivable that after a period of sev- 
i° yr, these objects will gradually develop a cometary 



eral 10' 

morphology with the youngest and most active star-forming 
regions located at one end of an elongated stellar body. 

5. Summary and conclusions 

Aiming at a better understanding of the different morpholog- 
ical subclasses of Blue Compact Dwarf Galaxies (BCDs) in 
the context of the evolution of star-forming dwarf galaxies, we 
have started a study of the "cometary" (iI,C ) subclass of BCDs. 
Contrary to the majority of BCDs where SF is confined to the 
inner part of an old circular or elliptical stellar LSB host, iI,C 
BCDs exhibit intense starburst activity close to one end of an 
elongated irregular LSB component. This intriguing morphol- 
ogy prompts two hypotheses: (i) iI,C BCDs are in fact dwarf ir- 
regulars (dls) observed during a major stochastic enhancement 
of their otherwise moderate SF activity and (ii) a set of phys- 
ical properties of the gaseous and stellar components favours 
the ignition of a burst with an amplitude comparable to that 
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typically observed in other BCDs at the outskirts of an older 
LSB component. 

For the two nearby cometary BCDs Markarian59 and 
Markarian 7 1 , a variety of deep ground-based and HST pho- 
tometric and spectrophotometric data are available. To inves- 
tigate the structural properties of the LSB component we de- 
rived surface brightness profiles (SBPs), corrected for system- 
atic effects which can be important at low S/N levels. Pho- 
tometric properties of the underlying host galaxy and the su- 
perposed starburst component were derived by fitting a sim- 
ple decomposition scheme to the SBPs; the radial extent of the 
starburst component as derived from this profile decomposition 
was verified using Ha- and colour maps. The spatial distribu- 
tions of heavy element abundances in the starburst regions and 
over the main stellar body were derived using long-slit spectra. 
After the superposed emission by ionized gas had been semi- 
empirically modelled and subtracted from the original spectra, 
a population synthesis analysis was carried out at several posi- 
tions along each slit to derive the properties of the underlying 
young and older stellar continuum. For Mkn 71 we also derived 
(B — V) colour-magnitude diagrams from HST data. The main 
findings of our analysis may be summarized as follows: 

1. The azimuthally averaged intensity distribution of the 
underlying LSB host galaxy of both iI,C BCDs can be ap- 
proximated by an exponential fitting law with a central surface 
brightness and a linear scale length intermediate between those 
typically inferred for iE/nE BCDs and dls/dEs. 

2. Spectral synthesis modeling of the starburst region and 
the main body implies the presence of an older population with 
a most probable age ^ 2 Gyr for the underlying host galaxies 
of Mkn 59 and Mkn 71. Ages up to 4 Gyr for Mkn 59 and 3 
Gyr for Mkn 71 are tenable within the model uncertainties, but 
higher values are unlikely. 

3. The average oxygen abundances were determined to be 
12+log(O/H)=8.0 (1/8 Z Q ) for Mkn 59 and 12+log(0/H)=7.8 
(1/14 Zq) for Mkn 71, which are typical among BCDs. In addi- 
tion, the metallicity distribution as derived for various elements 
in the vicinity of the starburst regions and along the major axis 
of the LSB-body shows only small scatter (~ 0.2 dex), suggest- 
ing that mixing of heavy elements has been efficient. 

The similarity of the two objects with respect to the ages 
and structural properties of their LSB components motivated a 
search for published data on other cometary BCDs. We found 
that five objects for which the required data are available share 
similar properties with Mkn 59 and Mkn 71 with respect to the 
results (1) - (2). These findings suggest that the specific star- 
burst morphology observed in iI,C BCDs comes along with dis- 
tinct physical properties of their LSB host galaxies, i.e. is not 
attributable to stochastic processes only. Hypotheses which we 
consider worth investigating are: 

(i) Cometary BCDs are relatively young objects, with ages not 
exceeding a few Gyr, thus systematically younger than "classi- 
cal" BCDs of type iE/nE. If true, the development of iI,C mor- 
phology may represent a late evolutionary stage of an iO BCD 
before it gradually assumes iE/nE-characteristics. 



(ii) The underlying host galaxies of iI,C BCDs with an age <; 1 
Gyr are moderately compact, in the sense that they show central 
surface brightnesses and exponential scale lengths intermediate 
between those typically derived for iE/nE BCDs and dls/dEs. 

The strongly extranuclear location of the starburst regions 
and the signatures of propagating star-forming activity along 
the main stellar body of an iI,C BCD suggest that the surface 
density distribution of cold H I gas in these systems resembles 
the optical morphology; published H I maps support this as- 
sumption. Spatially resolved interferometric studies will be of 
major importance for assessing the intrinsic processes regulat- 
ing SF in iI,C BCDs and exploring possible evolutionary links 
to "classical" iE/nE BCDs. 
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